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ABSTRACT 

Highly radiative expansion of a relativistic shell is shown to explain all observed features of 
the afterglows of the two bursts GRB 970228 and GRB 970508. In particular, in the first case 
the observed time-dependence oc of the soft X-ray flirx; is easily reproduced. The same 

model, when the surrounding matter density scales as a explains the afterglow of GRB 
970508, which may at first sight appear at odds with that of GRB 970228. In particular, it is 
shown that both the late peak in the optical luminosity and the flat time-dependence of the 
X-ray luminosity are simultaneously explained by nonuniformity of the surrounding matter, 
that the observed optical time-delay is correctly reproduced for standard parameter values, and 
that the time-delay and flux levels of the radio emission are also explained. 


Subject headings: gamma-rays: bursts - X-rays: transients - optical: transients 


1. Introduction 


It is an easy prediction of the fireball model that gamma ray bursts (GRBs) should show an afterglow, 
in the X-ray (Vietri 1997, Paper I) and at optical/radio wavelengths (Meszaros and Rees 1997). The energy 
source powering this luminosity is the kinetic energy of the (burn out) ejecta shell after the burst. This 
after-glow has become detectable thanks to an exceptional effort by the BeppoSAX science team which has 
managed to carry out TOOs just 8 hours after the burst GRB 970228, and then again about 48 hours later, 
thus providing the most complete coverage of the post-burst evolution of a GRB ever obtained in the X-ray. 

By combining observations of the Narrow Field Instruments (those of the TOOs) with both ASGA 
and Wide Field Gameras observations, Gosta et al. (1997a) were able to identify the main features of 
GRB 970228 as follows. First, the major burst lasts only « 5s; the afterglow begins just 30s after the 
major burst, has a time-dependence oc , with <5 = 1.32 ± 0.19, and the total fluence in the afterglow, as 
measured in just the 2 — 10 keV band, equals 40% of the whole burst (including the y-ray fluence!). 

I will show in the following that this behaviour can be easily understood within the model of Paper I. 
In particular, I shall show that the afterglow is due to the expansion of a highly radiative relativistic shell 
of matter, expanding in a constant density environment in the case of GRB 970228. 

A considerably different time behaviour is displayed by GRB 970508 (Gosta et al., 1997b) which 
shows an X-ray luminosity that decays very slowly (Gosta, talk delivered at the Elba Workshop on GRBs; 
Piro et al., 1997), and simultaneously an optical transient peaking well after the burst (Djorgovski et al.. 
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1997b). Qualitatively this behaviour is easy to understand. Since in the afterglow the dominant emission 
mechanism, synchrotron, is negligible for frequencies smaller than v^m the synchrotron turnover frequency, 
for the optical emission to set in one has to wait for relativistic effects to allow to move into the optical 
region. At the same time, the X-ray luminosity, which is well beyond is slow in decreasing. It thus 
appears as if the ejecta shell manages to remain relativistic longer than for expansion in a uniform medium 
{GRB 970228); expansion in a surrounding medium with decreasing density will be shown to provide a 
quantitative, as well as qualitative, explanation for the afterglow of GRB 970508. 

Understanding the behaviour of GRB 970508 seems especially important, in view of the minimum 
redshift determination of the optical transient (Metzger et al. 1997), which establishes the cosmological 
nature of GRBs: so long as ROSAT HRI observations linking precisely the site of the optical and X-ray 
transients are not forthcoming, even theoretical arguments linking the two may be valuable. 

In the next Section, I will discuss qualitative features in GRB 970228 which support the fireball model, 
while the hydrodynamics and X-ray emission will be discussed in Section 3. The optical and radio emissions 
are then discussed in Section 4. 


2. Qualitative deductions 

The light curve of GRB 970228 in both the 7 and X-ray bands is peculiar, in that it clearly shows 
a first peak lasting about w 5 s, a silence of « 30 s, and a new, longer peak lasting « 40 s. Costa et al. 
(1997a) attributed the second peak to the afterglow. There are three reasons for setting the second peak 
apart from the first one: first, the ■j/X ratio is radically different in the two events. Second, the luminosity 
of the second peak falls squarely on the extrapolation to early times of the time-dependence law that links 
the first and second TOOs to ASCA data. Third, the X and 7 light curves in the second peak start out 
nearly simultaneously, suggesting a common origin in a new shock, rather than the beginning of cooling 
in hot matter. Recognizing the second peak as part of the afterglow implies that the observed afterglow 
follows accurately a power-law time-dependence over nearly four orders of magnitude in time and five in 
flux, without apparently any band effect. 

The total soft X-ray (2 — 10 keV) luminosity released in the afterglow can then be integrated to show 
that it alone equals 40% of the whole burst luminosity (Costa et al. 1997a). Since the photon number 
counts go as cx; the total correction to obtain the bolometric luminosity must amount to a factor of 

a few, without changing the order of magnitude of the afterglow fluence. In short, the afterglow radiates 
about as much as the burst. 

This fact has two important consequences: first, the total energy radiated away by the expanding 
fireball is large. Thus, the expansion following the burst cannot be described by an adiabatic expansion, 
as was implicitly or explicitly assumed by other authors (Tavani 1997, Waxman 1997a, Wijers, Rees 
and Meszaros 1997), who treated radiative losses as a tiny perturbation to an otherwise adiabatic flow. 
Quantitative details pertaining to the radiative impulsive relativistic expansion will be derived in the next 
Section. 

Second, the approximate equality of fluences in the burst and in the afterglow makes it likely that the 
energy source powering the burst (the kinetic energy of the ejecta) is the same as that of the afterglow: 
if it were otherwise, i.e. if the the burst and afterglow were powered by different physical phenomena, an 
unlikely coincidence would result. 
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This equality is naturally accounted for in fireball models: when the Lorenz factor of the bulk expansion 
is low (ry w 100 ), the reverse shock is only marginally relativistic, and the total directed kinetic energy 
converted into internal energy, to be radiated as the burst, is 1/2 of the total energy budget (Sari and Piran 
1995). When, viceversa, ry is large, w 1000, the reverse shock is highly relativistic, and the ejcta are stopped 
dead in their tracks. In this case, the energy fraction converted into internal energy « 1, and the balance 
remaining as directed kinetic energy ^ 1, leading to a much weaker afterglow. This difference may account 
for the lack of afterglow from GRB 97011G or its relative weakness in GRB 970402. 


The following quantitative argument, supporting rj « 100, can be made. Of the two distinct mechanisms 
proposed to explain GRBs, for this low value of ?y only internal shocks (Meszaros and Rees 1994) can work, 
the other one, external shocks (Rees and Meszaros 1994) requiring much higher Lorenz factors rj rv 1000. 
In the interior case, two shells collide and stick at r ^ 10^^ cm, and from then on they continue a free 
expansion until they impinge upon the ISM, producing a second shock. This occurs at a radius 

1/3 . /1 nn\ 2/3 


^ sh 


= /- 
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ydTrnimpC^ry^ 


= 2.5 X 10^® cm eK^Ui 


( 1 ) 


where the total energy release is if = ifsilO®^ erg, and the ISM mass density is nimp cm~^. Thus, 
as seen from Earth the time-delay between the first and second shock and ensuing emissions is 
dt ~ Tsh/T.rf'c = 40 s(ry/100)“®/^. This value is essentially the time-delay between the first peak of the GRB 
970228, and the onset of the afterglow. Also, the simultaneity of the onset of X-ray and y-ray luminositites 
argues for a common origin in a new shock. I thus find the internal/external shock model for GRBs 
validated by the observed time-delay in GRB 970228, and the value of ry w 100 supported by observations. 


3. The X ray luminosity 


The time-dependence of the soft X-ray flux in a highly relativistic, radiative snowplow model was 
derived in paper I for late times (~ 1 month after the burst). I now derive the prediction for early 
post-burst times. 


The initial burst is contaminated by a baryon mass Mb = E/rjc^ with ry « 100; the shock with the ISM 
occurs at a radius given by Eq. |^. After this shock, the ejecta keep plowing through the ISM, shocking ISM 
matter and transforming its directed kinetic energy into internal energy. Assuming the post-shock cooling 
time to be short (to be checked later), the shocked matter piles up behind the shock in a cold shell, whose 
Lorenz factor 7 , for 7 3> 1, evolves according to (Blandford and McKee 1976) = —dM/M, where M 

is the total mass entrained by the shock, including the initial contamination. Assuming the surrounding 
matter density to vary with distance from the site of energy release r as p ex. r~°‘, it can easily be shown 


that 



and the time, as seen from Earth Ie, which scales as dE = dt/ 27 ^ : 

{7-2a)2p\ ((;^) 


drj^'^'^c, varies as 
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which, together, with Eq. || yields 
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where of course 


to = 


Tsh 


42 s 


100 \ 
V ) 


8/3 


(5) 


(7 - 2a)2r]'^c 7 - 2a 

The time to is the time-scale on which the scale-free solution sets in. Since we are mainly interested in 
times t ^ to, I shall from now on neglect the ±1 arguments inside the parentheses in Eqs. || and 

The parameter a can also be thought of as a ruse: for a = 3/2, we recover the results for the adiabatic 
shell. It should also be noticed that a is directly observable, as will be shown in the discussion following 

Eq. P. 


7-2q 9 X 10® S 7-2q 8 

tr = rj to = — --- 77 3 


The above solution is correct only for relativistic expansion; it breaks down for 7 Ri 1, ie., at a time tr 

7-2a ' 

which is very weakly dependent on a. The total time through which the shell remains relativistic is thus 
about a month, as shown in Paper I. 

The copious fluence in the afterglow phase (Costa et ai, 1997a) is evidence per se that the expansion 
must be highly radiative; however, this can also be checked through a simple, microphysical argument. 
For synchrotron emission, the synchrotron cooling time in the shell frame is ts = GTruiec/aTjaB^, to be 
compared with the local expansion time scale texp = r/"fc. I find 


TZ = 
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where I used an equipartition magnetic field B = ( 87 rnimpC^)^/^ 7 , and energy equipartition between 
electrons and protons, 7 e ~ mp'y/rrie. Using Eqs. |l| and I find 

/ y/(3-U-l 1 

K = '>^''(I) ^ ( 8 ) 


for a typical explosion, U 51 = 4 (Piran 1992). It can thus be seen that, for a = 2, 7^ <C 1, while, for a = 0, 
7^ = 1 for 7 = 2.8, close enough to the limit of validity (7 1) of Eq. to allow us to say that 77. < 1 is 

verified through the whole period in which the shell is relativistic. 

In the highly radiative limit, the total power radiated per unit time is given by a purely hydrodynamical 
argument (Blandford and McKee 1976). We have dE/dt = dTrr^cy^nimpC^, in the relativistic limit 7 1; 

in Earth time dE/dtE = 2'y‘^dE/dt. I find 


4q;-10 

7-2q 


Me 

The fraction of all power radiated in the X-ray region, fx, is computed in Paper I as 
fx = here the index p = 2.8 is the spectral index of the electron 

power-law energy distribution, dn ex 7 ^ 7 ^, with p fixed by the requirement that the observed spectrum, 

(X (Costa et ai, 1997a) be reproduced. Also, e„ and ei are the upper and lower limits of the BeppoSAX 
instruments, 2 and 10 keV respectively; em is the maximum electron energy. Proceeding as in Paper I, I 
find /x = (6 X 10“®/7)*^®“^^/^. In the end, I find 


Ex = 4.8 X 10 ® erg s ^ cm ^ 


ts 
1 s 


-s 
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The constant above is computed for a = 0, and a source distance D = 2 Gpc. It varies little with a, and 
given our ignorance of the source distance, it is not worth it to give its (complex) dependence on other 
parameters. Since 3 — p <C 1, the dependence upon waveband is very weak, explaining why SAX can observe 
such a strikingly extended power-law. 

The above discussion allows an immediate fit to the properties of GRB 970228. In this case, p = 2.8, 
and a = 0 (constant density) imply i5 = 1.38, in excellent agreement with the observed S = 1.32 ±0.19. This 
is clearly displayed in Fig. 1. GRB 970508 has a different time-behaviour, showing only a modest decline 
in flux by a factor 4.5 from O.Gd after the burst to 6 .Id after the burst (Piro et ai, 1997); this corresponds 
to an average decrease going like Fx oc Since the spectral properties of GRB 970228 and GRB 

970508 seem similar, we can use the equation above to see that, for a = 2, we have Fx oc in excellent 

agreement with observational data. It should be noticed that these data require a density gradient: for the 
highly radiative evolution, it was shown above that expansion in a constant density environment leads to 
Fx oc while adiabatic evolution leads to Fx oc where z w 0.9, the spectral slope of the burst, 

is most likely similar for GRB 970228 and GRB 970508. 


4. The optical and radio emission 


One may think at this point that the introduction of a declining density into which the fireball plows is 
but a simple minded way to account for the slowly declining X-ray flux, but I show here that this stratagem 
simultaneously explains why the optical luminosity is seen increasing, yielding quantitative predictions in 
excellent agreement with observations. 

Typical spectra of GRBs are of the form oc v°‘ for v < Vrm and oc for v > Vm, with a « 0 but very 
ill-determined, and 6 w — 1. Since the equipartition magnetic field B oc the turnover frequency as seen 
from the Earth I'm oc where jm oc 7, so that Vm oc The electron density in the shell frame Ue oc 7, 

and the comoving shell thickness dr = r/7, so that the comoving intensity oc /Bj^ oc rieBdr, 

and the observed flux as a function of observer time is oc oc (Wijers, Rees and Meszaros 

1997). From the above Eqs. ^and^I find F^^ oc F, where q — {2a — 3)/{7 — 2a). So long as Vm is shortward 
of the optical region, the optical luminosity is oc Fy^{uopt/i'm)‘^ oc F, where p = q + 4a(3 — a)/(3(7 — 2a)). 
For a ~ 0 and a = 2, I find p ra q^i 1/3, but ill-determined because of the uncertainty on a. 

Once Vm has entered the optical region, the luminosity must then decrease according to 

Fopt = ocF ; p={2a- 3)/(7 - 2a) ± 46 . (11) 

For a = 2, p = —1, for b = —1. The optical data from Table I are shown in Fig. 2, corrected for standard 
Galactic absorption and for a spectral shape oc i'~^. They show clearly that both the rise (less significantly) 
and the decline (more significantly) agree with the model for a = 2 . 

This model not only predicts that the optical luminosity should first increase and then decline, 
but it also reproduces correctly the time-delay between the beginning of the afterglow and the onset 
of optical emission. If the emission is pure synchrotron, there will be a minimum electron energy given 
approximately by 7 ^ = rnp'y/rrif, in the shell frame, emitting synchrotron photons at the turnover frequency 
= ^^m^B/AirmeC, with B given by the usual equipartition argument oc 7 . In the Earth frame, this 
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Table 1: Optical fluxes of GRB 970508. Time is measured from burst trigger, May 8.904 UT (Costa et al, 
1997b). _ 


logt 

(s) 

app. magn. 

m 

logv 

(Hz) 

band 

instrument 

ref. 

4.39 

«21.5 

14.74 

V 

Kitt Peak 

1 

5.04 

20.5 

14.74 

V 

Kitt Peak 

1 

5.04 

19.6 

14.56 

I 

Kitt Peak 

1 

4.41 

21.2 

14.64 

R 

Kitt Peak 

2 

5.00 

21.5 

14.92 

U 

La Palma 

2 

5.00 

21.0 

14.74 

V 

La Palma 

2 

5.00 

20.35 

14.64 

R 

La Palma 

2 

5.00 

20.2 

14.56 

I 

La Palma 

2,7 

5.26 

20.5 

14.92 

U 

La Palma 

2 

5.26 

20.3 

14.84 

B 

La Palma 

2,7 

5.26 

20.2 

14.74 

V 

La Palma 

2 

5.26 

20.1 

14.64 

R 

La Palma 

2 

5.26 

19.1 

14.56 

I 

La Palma 

2,7 

4.39 

21.33 

14.65 

Gunn-r 

Palomar 

3 

5.06 

20.17 

14.65 

Gunn-r 

Palomar 

3 

5.30 

20.15 

14.65 

Gunn-r 

Palomar 

3 

5.24 

19.65 

14.74 

V 

NOT 

4 

5.45 

20.53 

14.65 

Gunn-r 

Palomar 

5 

5.57 

20.76 

14.65 

Gunn-r 

Palomar 

6 

5.43 

20.9 

14.92 

U 

La Palma 

7 

5.43 

20.9 

14.84 

B 

La Palma 

7 

5.43 

20.6 

14.74 

V 

La Palma 

7 

5.43 

20.2 

14.64 

R 

La Palma 

7 

5.22 

19.8 

14.64 

R 

Loiano 

8 

5.46 

20.47 

14.64 

R 

M. Hopkins 

9 

6.33 

23.10 

14.64 

R 

HST 

10 


(1) Bond (1997); (2) Galama et al, 1997; (3) Djorgovski et al, 1997a; (4) Jaunsen et al., 1997; (5) Djorgovski et 
al., 1997b; (6) Djorgovski et at, 1997c; (7) Groot et al., 1997; (8) Mignoli et al., 1997; (9) Garcia et al., 1997; (10) 
Fruchter et al, 1997. 
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frequency is given by 
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47r 



7 

TTleC 


( 12 ) 


_1 /g 

As the shell decelerates, Vm will enter the optical region, « 2 eV. This occurs for ■^opt = ' , 

independent of all burst parameters except for the very weak dependence upon the density of the 
surrounding medium. From Eq. 7 = 7opt is reached at a time topt given by topt/to ~ 

Introducing numerical values, I find, for 0 = 2, topt = 3.4 days, in agreement with the observed topt = 2.5 
(Djorgovski et al., 1997b). In the case of GRB 970228 instead, corresponding to a = 0, I find topt = 3.9/i, 
well before the detection of the optical transient (van Paradijs et al., 1997), and thus in agreement with 
observations. 


Eq. 1^ is quite remarkable: not only does it show but a weak dependence upon the surrounding matter 
density and no other dependence upon burst luminosity, distance or beaming angle, it also tells us that, by 
mapping the times of flaring of the burst at different wavelengths, one can determine the time-dependence 
of the shell expansion, Eq. Thus the hydrodynamics of the problem at hand is, at least potentially, 
directly amenable to testing. 

It can easily be seen from Eq. that the turnover frequency cannot possibly have entered the radio 
region within the 5d time lapse within which the source was seen to flare (Erail et ai, 1997), so that it 
seems reasonable that this radio flaring is due to the source becoming optically thin. The optical depth 
below the turnover frequency for v < Vm, scales as r = Tm{v/vm'f' , due to the presence of a flat, most likely 
thermal, energy distribution of electrons (Tavani 1996). The optical depth at the turnover frequency 
can be computed approximately from the high frequency limit [u > Vm) as 


20 3 

3.1 X 10-® /tijX " 

- ^20/3 


(13) 


from which we find that the frequency which is becoming optically thin at any given time, r'ot = {jra)^^'^Vra 
is given by 

Vot = ^-^GHz[^^ . (14) 

The scaling has been chosen so to show that the detection of Erail et al. (1997), at 8.46 GHz a mere 5 d 
after the burst is correctly reproduced. Also, the corresponding expected flux level, Fot = , 

is given by 

4a-l 

/ f \ 3(7-2«) 

Fot = 0.6 mJy i—j , (15) 

fortuitously close to the observations of Frail et al. 1997, of 0.43 mJy. It should also be noticed that the 
inverted spectrum observed, oc which is so important in ruling out the hypothesis that the burst be 
due to a blazar, is characteristic of a radio source caught in the process of becoming optically thin. 


5. Discussion and summary 


An alternative and comprehensive model accounting for several properties of the afterglow of GRB 
970508 has been presented by Waxman (1997b), who postulates that the expansion is adiabatic. This 






makes the observed radiation small with respect to the afterglow energy, with the balance of this energy 
going into adiabatic losses; in fact, he postulates large {E w 10®^ erg) energy releases. There are two 
observational major differences between his model and the present one. First, in Waxman’s model, the slow 
decrease of the X-ray luminosity cannot be accounted for: in fact, Inverse Compton scattering (Waxman 
1997a) is incapable of decreasing the number of photons (and thus the X-ray luminosity) produced by 
synchrotron radiation, because it is well-known that the total optical depth to Thompson scattering in GRB 
ejecta shell is very small, w 10“® (Sari, Narayan and Piran 1996). Thus the model should show the same 
time-dependence of the X-ray luminosity as GRb 970228, which it does not (Piro et al, 1997). Second, 
since in his model the cooling time is long compared to the local dynamical time-scale, Waxman (1997b) 
expects an X-ray spectrum going as while I expect exactly like in the burst proper. Thus 

spectral observations will be able to tell whether the expansion is adiabatic (the X-ray spectrum is harder 
in the afterglow, Waxman 1997b), or highly radiative (afterglow and burst having similar spectrum, this 
paper). 

The main results of this paper are as follows: 

• The X-ray afterglow luminosities of both GRB 970228 and GRB 970508 are well-fitted by the 
deceleration of a radiative relativistic shell, plowing through external matter; in the first case, a 
constant density allows a good fit, while in the second one a power-like density distribution oc r“^, 
like that left over by a prior mass loss episode, is required; 

• the delay (« 30 s) between the first peak and the onset of the afterglow, in the X and 7 emission 
of GRB 970228 is quantitatively explained by assuming that the first peak results from an internal 
shock, and the afterglow from an external shock; 

• the existence of a delay between the optical luminosity maximum and the 7 peak in GRB 970508 
can most easily be explained by expansion in a nonuniform external medium, thus strengthening the 
interpretation of the X-ray light curve; 

• the optical time-delay for GRB 970508 is well reproduced as 3.4 days; 

• the late appearance of the radio flux, the observed flux level, and the peculiar inverted spectrum are 
all easily accounted for in this model. 


Thanks are due to Eli Waxman for showing me his work well before publication, to G. Ghisellini, G.C. 
Perola and E. Waxman for helpful scientific discussions, and to F. Pacini for organizing the Elba Workshop 
on Gamma Ray Bursts where this paper was completed. 
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Fig. 1.— 

X-ray flux in the afterglow of GRB 970228; data points are from Costa et ai, 1997a; the theoretical curve 
is from Eq. 

Fig. 2.— 

Optical observations of the afterglow of GRB 970508; data points are from Table I; straight lines (to be 
compared with the predictions of Eq. |ll|) are eye-fit to the data. 
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